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Abstract. We propose a novel scenario to generate primordial magnetic fields during inflation 
induced by an oscillating coupling of the electromagnetic field to the inflaton. This resonant 
mechanism has two key advantages over previous proposals. First of all, it generates a narrow 
band of magnetic fields at any required wavelength, thereby allaying the usual problem of a 
strongly blue spectrum and its associated backreaction. Secondly, it avoids the need for a 
strong coupling as the coupling is oscillating rather than growing or decaying exponentially. 
Despite these major advantages, we find that the backreaction is still far too large during 
inflation if the generated magnetic fields are required to have a strength of 0{W~^^ Gauss) 
today on observationally interesting scales. We provide a more general no-go argument, 
proving that this problem will apply to any model in which the magnetic fields are generated 
on subhorizon scales and freeze after horizon crossing. 
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1 Introduction 

Large scale magnetic fields of /uGauss strength have been observed in cosmic structures at 
a great variety of length scales and redshifts: from stars to galaxies to regions around high 
redshift quasars as well as from clusters and superclusters to the region of low density fil- 
aments [1-8]. Recently, different datasets have also been used to derive a lower bound of 
B > 10~^^ Gauss on the intensity of magnetic fields in the intergalactic medium [9-13]. 
There exist various astrophysical mechanisms of magnetogenesis and the amplification of 
weak seed fields by adiabatic compression during structure formation and by the galactic 
dynamo mechanism seem well understood [14-17]. Usually, primordial large scale magnetic 
fields are invoked to account for the seeds of nGauss strength needed for astrophysical am- 
plification. However, primordial magnetogenesis, active either during infiation or at phase 
transitions before recombination, is tightly constrained by Big-Bang Nucleosynthesis (BBN) 
and the Cosmic Microwave Background (CMB) [18-20] (in addition, such models do gener- 
ate non-gaussianities which can be used to constrain their model parameters [21-24]), and 
moreover most mechanisms have difficulties generating the required field strength or the co- 
herence length, see the reviews [25-31]. In fact, the required large correlation length is the 
most difficult property to explain. This applies equally well to astrophysical models and to 
primordial generation mechanisms. 

One promising focus for the genesis of large scale magnetic fields is cosmological in- 
fiation. Inflation is a quasi de Sitter stage of accelerated expansion of the Universe during 
which causally connected regions of space-time are stretched beyond the Hubble horizon. 
This enables the possibility of causal correlations on very large scales that re-enter the vis- 
ible horizon after the end of inflation, either during the radiation or matter eras. In order 
for any ampliflcation to occur during inflation it is necessary to step outside the Standard 
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Model (SM) comfort zone because of conformal triviality of electromagnetism (EM). Since 
EM is conformally invariant and de Sitter space (as well as Friedmann-Lemaitre-Robertson- 
Walker (FLRW) metrics in general) is conformally flat, this combination ensures that the 
Euler-Lagrange equations of motion (EOM) reduce to those in Minkowski space and no am- 
plification, or any other non-trivial dynamics, can take place. Therefore one has to construct 
extremely efficient couplings in order to overcome this adiabatic behaviour. In principle this 
can be envisaged but in most scenarios the created EM energy density gives rise to issues 
with backreaction on the inflaton dynamics. 

When a coupling of EM to the infiaton field (usually a scalar field) is implemented, it 
brings new terms in the EOM of the EM vector potential as well as the scalar field. For EM 
the coupling acts as an effective time-dependent mass term in the mode equation of the vector 
potential. However, as it grows the energy density of the EM field can grow accordingly, up 
to the point where it destabilises the de Sitter background. This can be shown to be a 
very general behaviour, and sets an upper limit on the maximal efficiency of a large class of 
mechanisms [32-35]. It turns out that in the most optimistic cases the final amplitude for 
the magnetic field, especially at large scales, is always several orders of magnitude below the 
desired level. Similar considerations apply to axial-type couplings [36] . The red thread which 
entwines these different mechanisms is the fact that all modes encompassed by inffation are 
amplified in a democratic way, from the lowest to the highest energies. The broadness of 
the amplified EM spectrum then makes it extremely hard to safely allow for very intense 
magnetic fields at some scales, as the total energy density has to be shared with all modes. 

In order to circumvent this issue, we explore a novel avenue in which magnetic fields 
are resonantly amplified during inffation thanks to an axial coupling to the inflaton ffeld^. 
Resonance happens in narrow bands such that only small ranges of modes can be amplified by 
this interaction. The resonance appears due to an oscillatory behaviour either of the coupling 
function itself or the dynamical evolution of the inflaton fleld^. The main purpose of this 
work is the analytic study of the consequences of such resonances and the requirements on the 
dynamics of the coupling. We demonstrate the feasibility of such a program by presenting a 
simple model at work. Despite the very different final EM spectrum, generation process and 
time, we will see that backreaction will still take place in a dangerous way before inffation 
can complete, and ultimately plague the ffnal fate of this class of models. 

We organise the presentation as follows. First, in Sec. 2 we make the case for the coupling 
between inffation and EM, outlining the reasons behind our choice (the Mathieu equation) 
and study the structure of the mode equation for the EM potential. Moving on to Sec. 2.2 we 
will both generally and explicitly treat the issue of building a working model. Sec. 3 is where 
we present the analytic estimates for the generated EM spectrum. The issue of backreaction 
is discussed in Sec. 4 and we wrap things up in the concluding Sec. 5. Appendix A lists the 
conventions we use for describing the quantised EM vector potential and its spectrum during 
inffation. Finally, appendix B contains a concise treatment of the Mathieu equation and its 
resonances. 

^Earlier work focused on the parametric resonance of gauge fields during the oscillatory final stages of 
inflation (preheating) and found insufficient amplification in typical scenarios, see [37]. 

^In general single field inflationary models with small periodic modulations in the potential [38, 39], an 
interesting consequence of the oscillating behaviour of the scalar field is the resonant enhancements of the 
correlation functions of cosmological perturbations, see for instance Axion Monodromy Inflation [40, 41]. It has 
been shown recently that such periodic modulations of the potential can lead to large resonant non-Gaussianity 
as well as other higher N-point functions [42-44]. 
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2 An oscillating coupling between the inflaton and electromagnetism 



2.1 Generalities and the Mathieu equation 

Our action is comprised of the free inflaton term £^ and the EM Lagrangian density Cem = 
— to which we add a chiral coupling to the inflaton 

S^EM = j d\^g ^-f{ip) FF , (2.1) 

where /{^p) is a dimensionless coupling function and F is the dual of the electromagnetic field 
tensor F, and FF = F^pF^^ = \ri"'^'^^ FapF^^ with rj'^l^'^^ being the totally antisymmetric 
Levi-Civita tensor with rj^^'^^ = -)-|(^|~i/2_ This coupling introduces an interaction term of 
the form fipFF in the EOM of the scalar field (we use the notation = df /dip). However, 
as we investigate the amplification of electromagnetic fluctuations we will neglect this term 
in a first approximation, returning to it when discussing the backreactions on the dynamics 
of inflation. 

We describe the electromagnetic field in terms of the vector potential in the Coulomb 
gauge (see appendix A for our conventions), (A^) = {0,Ai) with diAi = 0. The classical 
evolution equation of its Fourier transform in the helicity basis, A±{k, rf) reads 

^l^ + {k^ + hkf')Ah = G, (2.2) 

where h = ±1 stands for the positive and negative helicity, and /' = f,pip' is the conformal 
time derivative of the coupling function. This equation remains identical when we analyse 
the evolution of quantum fluctuations, see for instance [36]. Note that one may also consider 
couplings of the form /^((/j) F"^ or M{ip)'^A^ where /((/?) is a time-dependent coupling constant, 
and M{ip) is a dynamical photon mass, but we find it harder to resonantly generate magnetic 
fields in those cases. 

The central idea of this work is to employ a radically different coupling structure, which 
discriminates among different energy scales (and selects only a small portion of them), and 
is efficient in the perturbative regime of (2.1). These criteria can be met by an oscillating 
coupling term [c.f. Eq. (2.5)] which produces resonances. In the specific case at hand, Eq. (2.2) 
possesses resonant solutions as long as the coupling term is a periodic function of our time 
variable. We will focus on one specific case, the Mathieu equation, bearing in mind that any 
periodic coupling term would exhibit similar properties, and share qualitatively our results. 
Such an oscillating term with a constant, or slowly time-varying frequency, automatically 
selects and amplifies only a very small range of Fourier modes, leaving all other modes 
almost unperturbed. Furthermore, an EOM of the resonant type does not need any large 
parameter to produce appreciable enhancements. Our choice is also guided by our aim to 
be able to provide as much analytic understanding as possible, in such a way to clearly pin 
down the advantages, shortcomings, differences and similarities compared to the previous 
literature on the subject. 

Let us swiftly review the main properties of the Mathieu equation (more details are 
available in appendix B), which is typically written as 

d^y{z) 



dz2 



+ {p - 2q cos2z) y{z) = (2.3) 



Its solutions, the Mathieu functions, have well-studied instabilities in the (p, q) plane where 
they exhibit exponential growth. We can define the evolution variable as z = ur] for some 
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Figure 1. Left panel: The instability regions of the Mathieu functions are shaded in the {p, q) plane. 
The solid lines show where the stability properties change. The broken lines correspond to q — y/p 
(dashed) and q = \^ (dot-dashed). In terms of the dimensionless coupling constant, 7, of the generic 
oscillating coupling function given in Eq. (2.5) these represent the cases 7 = 1, 1/2. Right panel: The 
same chart but for our constrained parameters. Here p = (k/cu)'^ and q = jk/uj. 



constant frequency w which is hence a free parameter of the model; therefore we can obtain 
a Mathieu equation [c.f. Eq. (2.2)] by matching 

p = k^/uj^ , 2q cos 2ojrj = kf /uj'^ , (2.4) 

up to the sign of the coupling function which only determines the phase of the oscillating 
part of the solution. This means that both helicity modes are amplified almost identically, 
up to a phase shift. The choice of the Mathieu equation therefore automatically fixes the 
coupling function f{r]) to be of the form 

/(??)= 7 sin 2a;r/, (2.5) 

where 7 is a dimensionless coupling constant that cannot depend on k because /((/?) is only 
a function of the homogeneous background evolution, (p{r]). Before discussing the coupling 
as a function of if, let us mention that in our specific case we work with a constrained {p, q) 
set, since q = ^yk/uj = 7-^/^. In Fig. 1, left, we show the stability chart of the solutions 
to the Mathieu equation in the {p, q) plane. As examples, we draw the lines q = ^ and 
q = \ \/p through the chart that represent the restricted parameter spaces for the cases 7 = 1 
and 7 = 1/2. Also in Fig. 1, right, we use that restriction to draw the stability chart in the 
{k/oj,^) plane where one can clearly see that the resonance bands are around k = noj with 
integer n. 

If a mode A is resonantly amplified through the coupling to the inflaton, its amplitude 
roughly grows like |^| oc e^'^'', where ^ is the /c-dependent characteristic exponent or Floquet 
index, and oj is the frequency of the driving term in the Mathieu equation (see appendix B for 
details). For now this frequency will be kept constant. As we will elucidate later, this choice 
is physically not very natural and thus we will generalise to a time-dependent uj. Apart from 
the exponential growth the resonant solution also oscillates with constant amplitude. The 
effect of the opposite signs in front of the coupling functions for the two helicity modes of 
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the vector potential, A±, is a phase shift by 7r/2. As a consequence, the squared amphtudes 
of the two modes simply add up as 

+ ~ l^ipeSM-^^^, for A; -no;, n EN. (2.6) 

Here Ar] = rj — tji and Ai = A±{rii), assuming the same initial conditions for both helicities. 
The timespan Arj for which a given mode is resonantly amplified is therefore sensitive to the 
time when the process starts which would generically be the beginning of inflation. Notice 
that herein lies the unnaturalness of this choice: the final result depends critically on the 
initial time where we want to be able to set generic initial conditions. This will be alleviated 
when we allow for a time-dependent frequency which will also enlarge the range of modes 
that undergo resonant amplification. 

2.2 Background evolution 

Before we turn to a careful examination of the solution (2.6) and its consequences, we need 
to understand how a Mathieu equation can arise in the context of (successful) inflationary 
models. There are in general two ways to build a coupling function for which the dynamics 
of the inflaton result in the oscillatory behaviour that we need. 

Recall that we chose to stick to a coupling to EM of the form /(//) = 7 sin 2cj7/. The sim- 
plest way to obtain the correct time evolution is to simply solve for the background dynamics, 
which returns the functional form for 92(7?), and then engineer a coupling functional which 
translates this into the desired oscillations. This path leads to a f{(p) which automatically 
yields the Mathieu resonances through the change of variables 

An obvious disadvantage of this approach is that the functional form of f{^) may turn out to 
be quite unaesthetic. However, the point we want to make here is not one of attractiveness, 
but to show how in principle it is possible to build a working model of inflation which 
resonantly couples to EM. 

As an explicit example, we construct the coupling in the case of chaotic inflation with 
the potential V{ip) = ^m^ip'^ and the slow roll solution ip'^ ~ ArripN, where N is the number 
of e- folds until the end of inflation at rjf. A coupling of the form 

/((^) = 7 sin [2 CO r]f e('^/2™p)'-i/2l ^ (2.8) 

leads to the Mathieu equation and therefore the EM modes within the resonance bands are 
amplified exponentially as long as the slow roll approximation holds. 

Another solvable model is hybrid inflation with the potential V{ip) = Vo{l + g^V^^) 

where Vq sets the energy scale of inflation and m^ip^ <C for all values traversed by the 
inflaton. Inflation ends when the waterfall field becomes destabilised at (p = ipj, which may 
be chosen arbitrarily. Since ip{N) ~ (pj exp{N m'^ /iJip), we find the required coupling to be 

fip)=^sm[2u;^fiip/pfr"^/^"] . (2.9) 

The second possibility is to choose a simple coupling functional f{ip) and work out the 
inflaton Lagrangian implied by (p{ri) that is known from the chosen form of /(rj). Despite 
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the more attractive interaction term this simply transfers the comphcation to the form of 
dp. However, the inflaton field must roll down to the minimum of the potential and its field 
value must be a monotonic function of time. Similarly ip must always keep the same sign 
until reheating after the end of inflation. Hence both of these quantities can only contain 
a subdominant oscillating part, which is likely to make it difficult to construct an efficient 
coupling. On the other hand, the second derivative (p could oscillate so a coupling of the 
form /(□(/?) might work more efficiently, even though it would be necessarily suppressed by 
two more powers of the mass scale associated with the interaction. Due to these issues we 
focus on the first approach in this work. 

3 Resonant magnetic field generation 

With the background dynamics now specified, we are in a position to look at the evolution 
of the EM potential in some detail. In this section we first make analytic estimates of the 
resonant amplification in the simple case of a constant oscillation frequency. We then move 
on to the case of a time-dependent frequency which turns out to be more physical. 

3.1 Constant frequency 

If the frequency oo of the coupling is constant, only the Fourier modes within the resonance 
bands given by 7 and lo are amplified. In order to ensure that we are working in the pertur- 
bative regime of EM we take the coupling constant to be small, 7^1. Looking at Fig. 1 
we can see that for 7 <C 1 the resonance bands are very narrow in k/uj. The first band is 
located at fc = a; and is by far the most important and efficient. Its width can be estimated 
to be ~ 7 such that its boundaries are roughly at 

fc±^(l±|)^. (3.1) 

Within this first band the Floquet index can be estimated as 

^,u^[^^/^-{k/u-lf]"^ (3.2) 

such that we have ^k=uj — 7/2 at the peak. We only work with the first resonance band for 
the analytic estimates, as most of the EM energy density will be stored in this narrow band 
of modes. 

A resonant mode k € begins to grow exponentially according to (2.6) at the 

time inflation starts, r/j. The amplification lasts until k crosses the horizon at r]k = — 
because afterwards the cosine in the EOM does not go through any further period anymore 
and the oscillations of the coupling are effectively turned off. Hence, the time for which 
solution (2.6) describes the evolution of the mode k well is estimated as 

/^T] = r]k - rji . (3.3) 

Now we can estimate the power spectrum of magnetic fields amplified by the resonance. 
The power spectrum is oc k^ '^^\Ah\'^ , see appendix A, and the squared amplitudes of 
the resonant modes grow as exp{2fj,ujArj) according to Eq. (2.6). Since the resonance band 
is so narrow we can take the peak value for the Floquet index, ^ ~ 7/2, and model the 
fe-dependence by introducing a narrow dimensionless window function W{k/uj,'y) centred at 
k = oj that drops to beyond k±. Finally, we neglect % with respect to r]i in Ar/ which is a 
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good approximation for all modes that exit the horizon not too shortly after the beginning 
of inflation, as is the case even for modes which stretch up to 1 Gpc today. With these 
approximations we obtain the following spectrum at the end of inflation 



fc^|A(fc) P 

(27r)2a) 



Si = 4 W{k/u;,^) eM-l^Vfe^'f) • (3.4) 



after subtraction of the vacuum outside the amplified resonance band, Here a/ = 

a{r]f) and rji = rjje'^t with the total number of e-folds Nj and conformal time at the end of 
inflation rjf. The initial amplitude of the mode function, |A(A;)p, is given by the vacuum 
normalisation and has units 

The total EM energy density is given by f dln/c ((5^ + d^) where the electric field spec- 
trum is oc fe^l^^p. However, estimating the conformal time derivative of the resonant 
modes gives 

M^c^fiW\Ah\^^kMh\\ (3.5) 

and the contribution from the electric fields is negligible compared to that of the magnetic 
fields. For simplicity let us assume a standard Bunch-Davies vacuum initially, \Ai\'^ = {2k)~^, 
and replace the window function with a top-hat filter. Doing so we can integrate the spectrum 
to find the energy density at the end of inflation 



Pem^i[^ + ^) ^ exp(-7a;r,^ e^'f) . (3.6) 

However, we will see in the next section that backreaction will spoil inflation if this quantity 
is non-negligible. 

A very unattractive feature of this mechanism is the explicit dependence of the spectrum 
and the energy density on the initial time tji, and therefore the double-exponential dependence 
on the total number of e-folds Nj. The reason for this is of course the choice of a constant 
frequency co, which operates irrespective of the background dynamics. In addition to that, 
this frequency refers to periods in terms of the rj conformal clock, which is itself exponentially 
changing with N. In the next section we shall see how introducing a time-dependent frequency 
can mitigate this problem and give a physically more natural behaviour. 

3.2 Time-varying frequency 

We have seen that although a constant frequency does resonantly produce magnetic fields, 
the final result turns out to be sensitive to the total number of e-folds of inflation in a 
doubly exponential manner. In fact any direct sensitivity to the beginning of inflation feels 
unnatural and would bring about a dependence on the initial conditions of inflation, which 
is the opposite of what inflation is ordinarily invoked for. This may also lead to problems 
with the initial quantisation of the fields. The most obvious amelioration to the simplest case 
discussed above is to allow for a time-dependence in uj. Provided the time variation of the 
frequency is small compared to the oscillation time scale, the Mathieu equation will still be 
a good approximation to the resulting behaviour. 

In order to be able to set the initial conditions for all the modes in a consistent way, we 
choose the frequency 

u;{7]) = u e^"^ , (3.7) 
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Figure 2. The number of e-folds away from the end of inflation at which a mode k enters and leaves 
the resonance band, N{ri+{k)) solid and A^(ry_(fe)) dashed, are given for the model uj = ljc^"^. The blue 
curves correspond to uj = 0.2/Mpc, f3 = 10~''/Mpc; the purple curves to uj — 1/Mpc, (3 = 10~'*/Mpc; 
and the red curves are for oj = 5/Mpc, /3 = 10~'^/Mpc. We set 7 = 0.1 for all curves. The black 
dashed line is N{ri — —1/k) where the given mode exits the horizon scale. We zoom on the evolution 
of the mode k — 1/Mpc (solid black) for the second (purple curve) and third (red curve) choice of 
parameters: in the former case the mode becomes superhorizon before leaving the resonance band; 
the latter instead sees the resonance band fully traversed within the Hubble scale. 



where the strength of the time-variation is tuned with /3 while a) is the pivotal resonant 
frequency. In this case the resonance will not be operational throughout inflation for the 
same mode k, but will move across a range of modes. Thus, in addition to curing the initial 
condition problem, this allows for a more controlled amplification of a wider range of modes. 
This clearly makes it a more physical model. 

Given the simplicity of this model, we can make a number of reasonable approximations 
which will enable us to understand the resulting spectrum in an analytic way. The first step 
is to identify which modes resonate, and for how long. A mode k enters (+) and leaves (— ) 
the resonance band at r]± given by A; = (1 ±^/2)u}{r]±), c.f. Eq. (3.1). We find 



1 



In 



k 



(l±7/2)^ 



(3. 



We show these relevant times for three examples in Fig. 2. In this case a mode is resonating 
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for 



Ary = ry_-ry+ = -ln^^— ^J:.^, (3.9) 

where we used 7^1. Only for A; ~ a) do we need to worry about k exiting the horizon before 
leaving the resonance band. In this case Arj = r]k—rj^ ~ if^f^, where we can safely ignore rjk 
because typically f3 <^ lj. We observe that only the modes with (3 ^ k < u; pass through the 
resonance before horizon exit. Notice that the time spent in resonance band neither depends 
on k nor on the number of e-folds, a fact that reflects the more physical nature of the model 
with time-dependent frequency. 

With these estimates it is immediate to compute the amplification factor experienced by 
a resonant mode. The crudest approximation is obtained assuming that the Floquet index is 
a constant across the narrow resonance, ^ik — 7/2, and taking a; = A; at the time of resonance 
we find 

(the factor of 1 /2 in the exponent accounts for the fact that for A; ~ a) the resonance is active 
for only about half of (3.9)). Now that the coupling term in the mode equation is negligible 
well inside the horizon, we can safely set the Bunch-Davies initial conditions such that the 
spectrum is 5"^ oc A;^exp [7^/0/2/?]. A much better approximation can be found if we keep 
the full expression (3.2) for fi^ and integrate the exponent over time, which can still be done 
analytically, to find 

<5| = {r^'P - ^)^('^ -k)+e{k-u)^ , (3.11) 

where we have included the unamplified part for completeness, and the 9 are Heaviside unit 
step functions. This expression is quite accurate for u}/f3 < 10^ or so. The spectrum is very 
blue in the amplified region, k £ The two parameters of the model, /3 and uj, serve as 

natural IR and UV cutoffs in the theory. In order to have some amplification, /? should be 
small compared to the pivotal mode we want to amplify, uj. Let us finally compute the EM 
energy density by integrating the spectrum. 

pLi - exp [vr72(i/8/3] . (3.12) 

7r'^7^aj 

Similarly to the constant oj case we can safely neglect the electric field energy density and, 
therefore, the contribution to the total EM energy density comes only from the magnetic 
fields. 



4 Backreactions and subhorizon resonances 

The strong feedback from the generated EM energy density on the dynamics of the infiaton is 
a problem for many inflationary mechanisms of magnetogenesis. In our model there are two 
key properties that set it apart from most other scenarios: first, the resonant amplification is 
operational while a given mode is within the Hubble scale and freezes outside. Second, only 
a very narrow band of modes is amplified and not an entire power law spectrum. These two 
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characteristics compete for the final outcome of the backreaction calculation. A spectrum 
that is a narrow window at large scales reduces the total EM energy density as compared 
to that of a power law spectrum. On the other hand, the fact that the process ceases 
after horizon crossing demands a very large amplification factor in the first place. This is 
because the EM energy density dilutes like all the way from horizon exit until the end of 
inflation. We therefore need a huge magnetic field at horizon crossing in order for it to have 
an appreciable strength today. The larger the scale we consider, the earlier it exits and the 
more severe the problem becomes. 

Trying to balance these two effects, we find that backreaction always sets in too early 
and inflation is severely compromised. For instance, we can follow the growing magnetic 
field until it eventually catches up with the inflaton energy density. Or we can ask how long 
inflation could possibly last if we were to generate, at a given scale, a magnetic field to match 
the observed ones. Taking the second path, it is easy to have a rough idea of which parameter 
values would be needed in order to achieve a magnetic field that today peaks at 1 nGauss 
at 1 Mpc. Working with the time-dependent frequency (3.7), and fixing Co = 1/Mpc we find 
that the desired intensity 



We can now plug (4.1) into (3.12) to find that the EM energy density at the end of inflation 
is about 10~^/9(^ (where, for simplicty, p^p ~ V{ip)). Hence, it may appear that the generated 
magnetic fields are safely subdominant and inflation can proceed unharmed. However, in 
order to obtain these values at the end of inflation, we have to trace back the EM energy 
density to the point where the amplified window of modes left the horizon by taking into 
account the super horizon dilution of 1/a^; as we shall show in the next section, contrary to 
our expectations, the EM energy density greatly overtakes that of the inflaton at horizon 
crossing. 

4.1 No-Go for subhorizon amplification 

In conclusion, it is hard to generate cosmological magnetic fields at large scales strong enough 
to be measured today, or even to provide seeds for some late-time amplification mechanisms. 
This seems to be a manifestation of a more general pattern which is encountered in coping 
with subhorizon amplification mechanisms. Amplification of modes within (or around) the 
horizon and freezing outside appears to be a property of models with generic axial couplings 
even without resonances, the reason being that the coupling is suppressed by a factor k 
outside the horizon, see [36]. 

We will present a more formal "proof" of this fact here. The assumptions, encompassing 
the models we have analysed so far, are: 

1. the amplification mechanism is only active inside the horizon, and ceases, with conse- 
quent freezing of the vector potential, once a mode crosses the horizon; 

2. the only effect of the reheating stage following inflation is to alter the overall decay of 
the total inflaton energy density and the evolution of the scale factor, and has no direct 
impact on the EM field. 

Let us start with the aim of producing a magnetic field strength i?o* today at a correlation 
scale k-t . To minimise backreaction we minimise the total magnetic energy density by consid- 
ering a power spectrum that is a Dirac delta peak, Pg = Ak^,5xi{k/k^ ~ !)• Then the energy 



Or 




g^7''i/i6/3 ^ nGauss ^ 1.75 x 10"^^ Mpc 



-1 



(4.1 
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density is = A/c^/(27r)^ and the coefficient A is chosen such that = ^-Bq*- ^^^^ 
neglect the energy density in the electric field during inflation, in order to be as conservative 
as possible. We can then calculate the energy density of the magnetic fleld N e- folds before 
the end of inflation 

„ ( Ar\ — ^endAN _ n -4 AN (a r,\ 

Pb[JS)-Pb e - p^a^^^e , (4.2 j 

where aend is the scale factor at the end of inflation (normalised today) and pIj"'^ is the 
magnetic energy density at the end of inflation. We have assumed that the behaviour is 
obeyed throughout the evolution of the post-inflationary Universe as well.^ 

The backreaction problem is caused by the fact that pB decays much more quickly 
than the energy density of the inflaton. The simplest and most optimistic way to treat the 
evolution of the Universe after inflation is to assume it was always radiation dominated. 
In that case the ratio fi^ = pb/ Ptota.i is constant after inflation ends. For a conservative 
magnetic fleld strength of i?o* = 10^^^ Gauss today we have ~ 10^^^ <^ 1, and for a 
more realistic calculation it would clearly be larger. However even for this optimistic case, 
we can see that, evolving the magnetic field backwards into infiation, its energy density will 
dominate after only 

iV* = -^InOs f« 13 (4.3) 

e-folds if we assume a de-Sitter expansion, independent of the scale k^, and the inflationary 
energy scale. More realistically, in slow roll the inflaton energy density slightly increases 
going backwards into inflation which improves the situation marginally. However, the obser- 
vational bound on the tensor-to-scalar ratio means that the energy density cannot increase 
very significantly. This clearly demonstrates that we could only have magnetic fields with 
an interesting amplitude (i^o* ^ 10^^^ Gauss) on very small scales, k > = ^{N^), cor- 
responding to scales which crossed the horizon as late as < 13 e-folds before the end of 
infiation. The problem is alleviated a little with a lower energy scale for infiation, as large 
scales would exit the horizon at smaller A'^, but the improvement is minimal. 

One may turn this reasoning around, and ask what is the maximum magnetic field 
strength which could be produced within the limits set by the backreaction. In this case we 
can fix the ratio Q.b{N) = 1 at the time of horizon exit, which is the earliest time at which 
we have absolute confidence in the time-evolution of the magnetic field (without specifying a 
concrete model). From then on and until the end of inflation this ratio dilutes away as e~^^: 
the flnal ratio, for large scale modes which typically cross the Hubble scale around > 50 
is therefore Qb ^ 10"^'^, and the corresponding magnetic field today would be Bq <^ 10"^'' 
Gauss. Hence, any reasonable duration for inflation automatically leaves us with only crumbs 
of magnetic energy today. 

Finally, we should emphasise that this no-go theorem is not specific to the case of a 
narrowly peaked magnetic power spectrum. As mentioned above, assuming a delta function 
for the spectrum simply is the most conservative choice and any broader distribution or power 
law spectrum gives even more stringent constraints. 

^This assumption is well justified since we are looking at very large scales, and there is no risk of incurring 
any (direct or inverse) cascades. Even if an inverse cascade was active it is unlikely that it could have a huge 
effect, see e.g. [36]. 



- 11 - 



5 Summary and conclusions 

In this work, we have made a novel proposal for the generation of primordial magnetic fields 
induced by an oscillating axial interaction between the inflaton and the EM field. This may 
lead to a highly efficient resonant production of magnetic fields during inflation, and, unlike 
other scenarios, without the need for a coupling whose value grows or decays dramatically 
during inflation. Through studying the Mathieu equation, we can see that only certain 
wavelengths, which correspond to instability bands leading to exponential growth, will be 
amplified, which is another key advantage of our scenario. While most alternative scenarios 
ultimately lead to a blue magnetic field spectrum because of the onset of the redder and 
stronger electric field, in our scenario one may tune the parameters of the coupling so that 
only modes of observational interest are amplified, leaving the electric field modes behind; 
this leads to a sharply peaked spectrum, which will spread over a more realistic range of 
wavelengths by the further evolution in the plasma. 

We have proven how for any infiationary potential one may construct a coupling such 
that the magnetic fields experience the resonances of the Mathieu equation, both through 
a general formula for reconstructing the coupling and through two explicit examples. The 
results are rather complicated. However, we only chose to work with the exact Mathieu 
to be able to provide analytic estimates of the generated fields and to neatly capture the 
underlying physics. Notice that any periodic coupling would give rise to resonances and 
similar spectra and we therefore do not have to restrict to exactly those couplings derived 
in Sec. 2.2. Alternatively, and even in the case of the Mathieu equation, one may choose 
a simple coupling and instead engineer the dynamics of the infiaton to generate the same 
resonance. 

We have studied both the case of a constant resonant frequency, which is analytically 
simplest but sensitive to the beginning of infiation, and the more realistic case of a slowly 
varying coupling. In the case of a constant frequency, the modes in a resonance band will 
grow exponentially from the beginning of inflation until horizon crossing, leading to a doubly 
exponential sensitivity (in terms of the number of e-folds) on the initial time of inflation, and 
to theoretical embarrassment in flnding a consistent set of initial conditions. We have shown 
how a slowly varying frequency solves these problems. 

Generating primordial magnetic flelds with large correlation scale and without backre- 
acting is notoriously difficult. We have outlined above several reasons why our new method 
should offer significant advantages. Unfortunately the backreaction does remain a severe 
problem in our scenario, even when the magnetic field is concentrated in a very narrow band 
and negligible everywhere else. The issue is that the magnetic field generation only takes 
place while the chosen modes are still subhorizon, and after horizon crossing they freeze. 
Hence, their energy density is diluted like after crossing the horizon. This is not so 
severe after infiation until today, since the radiation and matter densities also decay rapidly. 
However, the effect during inflation is dramatic because the inflaton energy density is nearly 
constant by construction, and so the relative energy density of the magnetic field will be 
diluted exponentially in terms of the number of e-folds from horizon crossing until the end 
of infiation. 

Under the most optimistic assumptions, we have traced a single Fourier mode of strength 
10~^^ Gauss today backwards into inflation. We found that already = 13 e-folds before 
the end of inflation backreaction sets in, compromising the course of slow roll inflation. It 
is therefore clear that our mechanism, despite its advantages, cannot be invoked to generate 
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large scale magnetic fields of sizeable strength today. 

We emphasise that this no-go result does not only hold for the mechanism studied here. 
It can directly be applied to all mechanism of inflationary magnetogenesis where the vector 
potential modes are constant outside the horizon. Only special reheating scenarios or inverse 
cascades in the plasma era that change the shape and/or amplitude of the magnetic spectrum 
in an extreme way could potentially mitigate this result. In conclusion, the problem of the 
origin of cosmic magnetic fields again remains wide open. 

Note added. After this work appeared on the arxiv, a quite general no-go result for 
primordial magnetogenesis was presented in [45, 46]; in which it was shown that the magnetic 
fields generated during inflation typically yield unacceptably large anisotropic stresses during 
the subsequent radiation era. It has not been checked whether this result will apply to the 
model proposed in this paper. 
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A Conventions 

We work with a metric signature ( h + +). For tensor components, Greek indices take 

values ... 3, while Latin indices run from 1 to 3. We employ Heaviside-Lorentz units such 
that c = h = ks = €o = fiQ = 1. The reduced Planck mass is defined as nip = (SvrG)"^/^. 
We normalise the cosmic scale factor to unity today so that comoving scales become physical 
scales today. 

For the EM fields, our conventions follow [47]: 

= F^^u'' , = K^^^p^F^'Pu^ = F^au'' , (A.l) 

where Eau'^ = = B^u'^. In a perturbed FLRW metric an observer has the four- velocity 
ti^ = a^^(l, 0) + 0(1). As a consequence one finds (in Coulomb gauge) 

{E^) = (o, -^A^) , {B^) = (o, KjkdjAk) . (A.2) 
Introducing the orthonormal spatial basis 

with \ef\^ = 1, k = k/k, (A.3) 

and the helicity basis 

ei = (e? ± zsl) , (A.4) 

the vector potential takes the form +A-e-. After quantisation of the vector potential, 

we can study the evolution of the Fourier modes, Ah{r], k), with respect to the helicity basis 
for the polarisation states, /i = it. 
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If the magnetic field generated by some process is statistically homogeneous and isotropic, 
its spectrum is determined by two scalar functions Psik) and PA{k). Since the magnetic field 
is a divergence- free vector field the two-point function of the Fourier components of the mag- 
netic field can be written, for a comoving wave vector k, as 



{Bi{7],k)B*{r],q)) = ^^6{k - q){{8ij - k;kj)Ps{v, k) - le^jnknPAiv, k)} , (A.5) 

where Ps and Pa are the symmetric and anti-symmetric parts of the power spectrum, re- 
spectively. The symmetric part of the spectrum determines the energy density: 

{Biirj,k)B*{ri,q)) = {27rf 6 {k - q)Ps{ri, k) . (A.6) 

With respect to the helicity basis the spectra can be directly written as 

Ps/Airi,k) = kmA+{v,k)\^±\A^iv,k)\') , (A.7) 

where the upper sign corresponds to Ps and the lower sign to Pa- Here we use the non-trivial 
result widely applied in inflationary cosmology that at late times, the vacuum expectation 
values of the fields generated during inflation can be interpreted as stochastic power spectra. 
Finally, sometimes in the literature the power spectrum is identified with 

~^B^J^2Psiv,k) ^ ^B = j^,Ps{ri,k), (A.8) 

which is often denoted as dps /dink = 5'^. When integrated 

dk 

^61, (A.9) 
it gives the average magnetic energy density, pB = \B^ . 

B The Mathieu equation, stability and characteristic exponents 

The Mathieu equation we are working with is written in the form (2.3): 

d^y[z) 



dz2 



+ {p - 2qcos2z) y{z) = {). (B.l) 



Stability chart: In order to find in which regions of the {p,q) plane (recall that p is 
essentially a normalised k"^) the solutions to the differential equation diverge, we need to 
identify the solutions with period vr or 27r. It turns out that such solutions are always 
expressible in series of sines or cosines such as 

Period vr / = E^=o «« cos 2nz 
Period 27r / = E"=o cos(2n + l)z 

Period 2vr | = ^ gi^(2n + 1). ' ^^'^^ 
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Plugging these back into the original Mathieu equation, we find a system of recursive relations 
such as 



System A 

System B 
System C 
System D 



= (n > 2) 
(n > 2) ' 



(B.4) 

(B.5) 
(B.6) 
(B.7) 



pao — qai = 
{p - 4)ai - q{2ao + 02) = 
{p - 4n2)a„ - q{an-i + a„+i) 

{p - 4)61 -qb2 = 
{p - 4:V?')bn - q{hn~i + hn+i) 
(1 -p + q)co + qci = 

[{2n + 1)2 -p\cn + q{Cn-l + Cn+l) =0 (n > 1) 
{l-p- q)do + qdi=0 

[{2n + 1)2 -p]dn + q{dn-i + d„+i) =0 (n > 1) 

where we have used a few handy multi-angle trigonometric relations such as 

2 cos 2z cos 2nz = cos 2{n + l)z + cos 2(n — \)z (n > 1) , 
2 cos 22 sin 2n2; = sin 2(n + l)z + sin 2(n — l)z (n > 1) , 
2 cos 2z cos(2n + \)z = cos(2n + 3)2; + cos(2n — \)z (n > 0) , 
2 cos 2z sin(2n + l)z = sin(2n + 3)z + sin(2n — l)z (n > 0) . 

The systems (B.4), (B.5), (B.6), and (B.7), must have zero determinants in order to 
possess non-trivial solutions for at least some of their coefficients a„, 6„, c„, and dn- Therefore 
the equations 





Da 



Db 



Dc 



Dd 



p -q 

-2q p — A —q 
—qp — lQ' 

p-4: -q 
—q p — 16 —q 
—q p — 36 



1 — p + q q 
q 9 — p q 
q 25-p 

1 — p — q q 
q 9 — p q 
q 25-p 



0, 



(B. 



0, 



(B.9) 



0, 



(B.IO) 



0, 



(B.ll) 



must be satisfied (not simultaneously) . Each time one of those is satisfied we have identified 
a solution with the required period, that is, a boundary between stable and unstable regions 
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P 

Figure 3. Stability chart: the bands alternate between stable and resonant regions, from left to right 
or vice-versa, beginning with a stable one on either side. 

in the {p, q) plane. We will confine ourselves to the case where {p, q) G M. The stability chart 
is depicted in figure 3. 

Characteristic exponents: We now seek some more solutions in the form (Floquet's) 

F^z) = e^^'Piz) , F,{-z) = e-^^'Pi-z) , (B.12) 

where P is a regular oscillating function of period vr. Thanks to Floquet's theorem, as long as 
v ^ Z the two solutions F{z) and F{—z) are linearly independent, and the general solution 
to (B.l) can be written as 

y{z) = AF,{z) + BF,{-z), (B.13) 
with A and B constants. Again, as long as u ^Ij we can decompose Floquet's solutions as 

+ 00 +00 

Fu{z)= /ne^^'^+'"^'andF,(-z)= /ne-'("+'")^ (B.14) 

n=— oo ?i=— oo 

Obviously, the region of instability we are interested in is that for which v = i^i with /i real. 
Using the same machinery as before we could in principle find the solutions for /i on the (p, q) 
plane. In practice, we plug the expansions (B.14) back in the differential equation, extract 
recursive equations among coefficients, and their corresponding determinants must vanish. 
The resulting master recursive equation is the same for both solutions, and reads 

+00 

Y { [(^ - 2m)2 +p\ U - q{f^_, + e^'^^i^.n. ^ ^ (B -^5) 

n=— oo 

from which the recursive relation 

[il^ - 2inf +p]fn- qifn-l + fn+l) = , (B.16) 
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can be extracted. In principle one then analyses the determinant as before, reading 



where we define C = <?/ [(^/^ + 
summation 



■P 



1 C 
c 1 

and N - 



(B.17) 



AfxiV 

oo. This determinant can be recast as the 



D 



(JV) 



N 



a, 



.N 



a 



,N-1 
I 



a 



.N-2 
1-2 ' 



1 , a 



N 



0. 



(B.18) 



1=0 



As far as we are aware of, the recursive relation can not be solved: since we need Dp = 
limAT^oo -Cp^^ then we are not able to obtain an exact expression for the determinant. 

There is another path one can follow in order to have exact expressions which one can 
use to plot the characteristic exponent. Let us take two solutions of Mathieu's equation 
which we call yi and y2 such that 



yi(o) - 
y'i(o) ^ 

2/2(0) 



1 



: 
: 1 



(B.19) 

(B.20) 



Now, the following relations can be proven: 



arccos [yi (tt)] , 

arccos [l + 2y;(7r/2)y2(vr/2)] • 



(B.21) 
(B.22) 



Since we ask that iJ, is real, we are only looking at purely imaginary values for the arccosine. 
In turn this means we are only interested in its argument being larger than +1 (or smaller 
than -1, in which case the real part is just vr and simply causes a sign flip in Fi,{z)). The full 
explicit forms for yi and y2 are combinations of sine and cosine Mathieu functions (se and 
ce, respectively): 



while 



y2{z) 



ce'{a, q, O)se(o, q, z) — se'{a, q, 0)ce(a, q, z) 
ce'{a, q, 0)se{a, q, 0) — se'{a, q, 0)ce(a, q, 0) ' 



ce(a, q, z)se{a, q, 0) — se{a, q, z)ce{a, q, 0) 
ce'{a, q, 0)se(a, q, 0) — se'{a, q, 0)ce(a, q, 0) 



(B.23) 



(B.24) 



After choosing a value for q, it is simple to plot the characteristic exponent as a function 
of p. In the context of the EM mode equation with an axial coupling to the inflaton these 
parameters correspond to p = k'^ /uj'^ and q = jk/u}. Figure 4 shows examples of the resonance 
bands for this case. 
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Resonance Bands with 7 = 1 
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Resonance Bands with y = 1/2 
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Figure 4. Characteristic exponent and resonance bands. The vahie of q is y/p (left panel) or y/p/2 
(right panel), that is, 7 = 1 (left) or 7 = 1/2 (right). 
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